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Abstract. - The Ising-like spin ice model, with a macroscopically degenerate ground state, 
has been shown to be approximated by several real materials. Here we investigate a model 
related to spin ice, in which the Ising spins are replaced by classical Heisenberg spins. These 
populate a cubic pyrochlore lattice and are coupled to nearest neighbours by a ferromagnetic 
exchange term J and to the local (1, 1, 1} axes by a single-ion anisotropy term D. The near 
neighbour spin ice model corresponds to the case D/J — » oo. For finite D/J we find that the 
macroscopic degeneracy of spin ice is broken and the ground state is magnetically ordered into 
a four-sublattice structure. The transition to this state is first-order for D/J > 5 and second- 
order for D/J < 5 with the two regions separated by a tricritical point. We investigate the 
magnetic phase diagram with an applied field along [1, 0, 0] and show that it can be considered 
analogous to that of a ferroelectric. 



Background. - Geometrical frustration represents a recipe by which condensed matter 
can be disordered even in the absence of substitutional disorder jl| . The canonical example 
is the proton disorder in ice, which was famously shown by Pauling to be the origin of the 
experimentally observed ground state entropy |§,p|. Anderson later illustrated a direct map- 
ping of Pauling's model onto the Ising antiferromagnet on the pyrochlore lattice Q. However, 
Anderson's antiferromagnet does not appear to occur in nature and the mapping was for 
many years something of an academic curiosity, albeit one that inspired interest in frustrated 
magnetism ^Q). Recently, however, this has changed with the presentation of an alternative 
mapping named "spin ice", motivated by experiments on the cubic pyrochlore H02T12O7 S. 
This material is closely approximated by a pyrochlore ferromagnet with (1,1,1) Ising-like 
spins, which maps onto the ice model || and hence onto Anderson's antiferromagnet P,pl|. 
The behaviour of Ho2Ti 2 07 || and the related compound Dy 2 Ti 2 07 |l^,[l4| is described to 
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Fig. 1 - The spin ice mapping a) Coordination of H2O molecules in ice represented by proton 

displacement vectors, b) Spins in the spin ice model. 

a good first approximation by the spin ice model and has recently been shown to be repro- 
duced to an excellent degree by a related "dipolar spin ice model" including long range dipolc 
interactions [pjl|l6l[l8|]. These materials are part of a series of rare earth titanates that are in 



In Ho2Ti207 and Dy2Ti2C>7 




general of interest for their frustrated magnetic behaviour 

The mapping may be illustrated as follows (|] (see Fig. 
the magnetic rare earth ions occupy a cubic pyrochlore lattice, an array of corner-linked tetra- 
hedra. This lattice is the medial lattice (i.e. the lattice formed by the mid-points of the 
bonds) of the diamond- type oxide sublattice in cubic ice Q| . The proton disorder of ice can be 
described by displacement vectors that populate the vertices of this lattice. The "ice rules" , 
the condition that there are two protons near to, and two further away from each oxide, cor- 
responds to the condition that two vectors point into and two point out of each tetrahedron. 
In Ho2Ti2C>7 and Dy2Ti2C>7 the rare earth ground state is an effective S = 1/2 doublet with 
local (1,1,1) quantization axis and so can be represented by a classical Ising-like vector anal- 
ogous to a proton displacement vector. Net ferromagnetic coupling between near neighbour 
spins completes the analogy with the ice model, corresponding to the frustrated proton-proton 
repulsion imposed by the ice structure. The spin ice model describes qualitatively the results 
of many experiments on these materials, including the absence of magnetic order down to the 
lowest temperatures [SClLO. 

The Model. - The macroscopic degeneracy of the spin ice model makes it very susceptible 
to perturbations. That is, one can expect correction terms to the idealized Hamiltonian to 
dramatically change the low temperature properties by lifting the degeneracy, moving the 
system into a long range ordered state. This is the case for dipolar interactions, as long as 
thermodynamic equilibrium is maintained. In this paper we pursue the effects of perturbation 
by studying a continuous spin ice model, with classical Heisenberg spins replacing the discrete 
Ising-like degrees of freedom. Our model is defined by the spin Hamiltonian 



H = -J S l .S J -Dj2(S t -d i ) 2 ~Y. li A 



(1) 



<i,3> 



where the Si are classical vectors of unit length. The cubic crystal fields of magnitude D 
are along the four (1, 1, 1) directions for the four spins of the primitive unit cell. The extra 
refinement of continuous, rather than Ising spins would be expected to capture the effect 
of quantum fluctuations, which must be present in some extent in all real materials. One 
would also expect it to be an accurate starting Hamiltonian for any prospective ice or spin 
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Fig. 2 - a) Monte Carlo simulation of the magnetization versus temperature for Hamiltonian |l| 
Thermal hysteresis was observed for large D/J; the results shown are cooling curves at (from left 
to right) D/J = 22.91, 15.49, 10.47, 7.08, 4.79, 3.24, 2.18, 1.0. b) Spin structure of the q = basis for 
D/J — 7, illustrating the canting of the spins towards [0, 0, 1]. 
direction. 



ice like materials with an exchange constant J in excess of the 0(1°) K temperature range 
of Ho2Ti20y and Dy^T^Oy, where dipolar interactions are important. The Hamiltonian (^) 
describes the transition from spin ice {D/J — ► oo) to ordinary ferromagnetic (D/J — ► 0) 
behaviour. We have investigated it by Monte Carlo simulation using the standard Metropolis 
spin- flipping algorithm for systems of size N — 432 — 11634 spins. All simulation lengths were 
100,000 Monte Carlo steps per spin with 30-40,000 steps used for initial equilibration. 

Spontaneous Magnetization. Our first result is illustrated in Fig. 0, where we plot 



Monte Carlo simulation results of the scalar magnetization per spin m = 




as a function of temperature T/ J for < D/J < 25. For all the finite D/J investigated 
there is a transition to an ordered state with q — [0,0,0] propagation vector. It is clear 
that the transition only disappears in the spin ice limit, D/J —> oo and therefore that the 
infinite ground state degeneracy is lifted, for all finite D/J, in favour of a long range ordered 
q = state. The transition is marked by a maximum in specific heat and susceptibility. It 
is second-order for D = (the simple ferromagnetic case) but clearly first-order for large 
D/J, suggesting the presence of a tricritical point [|l9],|2(J at intermediate D/J. To locate 
the tricritical point Dt c / J we have performed a finite size scaling analysis of the maximum in 
the susceptibility Xmax- Xmax should scale as L 1 ^ in the second-order regime ^l), where 7 
and v are the critical exponents describing the divergence of the susceptibility and correlation 
length. Assuming the Ising universality class, one has 7/f = 1.969 ~ 2 p2fl . Such L 2 scaling 
was observed for D/J < 5 only, and a careful analysis suggests that the tricritical point occurs 
at D tc /J = 5.0 ±0.5. 

To further understand the results of Fig. [|, consider the nature of the ordered state: a 
q = state has a magnetic unit cell consisting of a single tetrahedron. For the Ising system 
(D/J — » 00) the spins are directed along the local (1,1,1) shown in Fig. [j]b. We 

call this the "(1,1,1)" state which has energy per spin Emii) = —J/3 — D and overall 
magnetization m = 1/V3 along z = [0,0,1] . The state is frustrated by the crystal field 
with, for each tetrahedron, four bonds making a contribution of —J/3 and two making a 
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contribution of J/3. Introduction of continuous spins allows the system to regain frustrated 
bond energy at the expense of field energy by relaxing away from the body centered cubic 
directions; see Fig. ||b. In disordered spin ice ground states this relaxation is itself frustrated 
and it is easy to see that the maximum energy gain will occur for relaxation away from the 
ordered q = state. For finite D/J we therefore expect degeneracy to be lifted in favour of 
this state and to have a canting of the spins towards the [0,0, 1]. In the limit D/J — > one 
should obtain a perfectly aligned ferromagnetic (FM) state, which we call the "[0, 0, 1] state". 
For finite D it has energy -E[o,o,i] = — 3 J — D/3. 

We have estimated the ground state spin orientation for large and small crystal field 
by making the reasonable ansatz that the relaxation away from, or towards the (1,1,1) is 
homogeneous for the four spins of the unit cell. Comparing with numerical data this seems a 
posteriori to be correct. For small a = D/J a small relaxation away from [0, 0, 1] and towards 
(1,1,1) gives for the magnetization m and the energy E per spin 



This result illustrates the fact that the energy and magnetization change very slowly from the 
ferromagnetic value for D increasing from zero to D ~ 0{J). In the opposite limit of strong 
crystal field, we define 7 = J/D and allow a small relaxation away from the (1,1,1) state 
towards the [0, 0, 1] state. In this case we find for small 7 



which again suggests a slow change in magnetization with D/J. 

To test these expressions, the magnetization was computed by performing zero temperature 
Monte Carlo simulations. In Fig. || we plot these estimates versus the anisotropy D/J 
and compare with the small anisotropy result (Q) and the large anisotropy result (|3|). The 
agreement is seen to be excellent in both cases, with each asymptotic equation breaking down 
near the tricritical point at Dt c / J = 5. It is therefore tempting to associate the change in 
order of the transition with a non-linear cross over from a regime that is nearly ferromagnetic, 
to one that is nearly q = spin ice. 

Magnetic Phase Diagram. - Further insight into the nature of this cross over is gained by 
a consideration of the magnetic phase diagram. We consider the behaviour of the continuous 
spin ice model in an external magnetic field applied along the [1, 0, 0] direction which sustains 
the symmetry of the ordered state for all finite D/J. In ref . ( |l3[ ) we showed that the magnetic 
phase diagram of the near neighbour spin ice model with the field along [1, 0, 0] has two lines 
of first-order phase transitions that separate phases of different magnetization. Both lines 
terminate in critical points, and we have argued that the phase diagram can be considered 
analagous to that of a liquid-gas system, where a line of first-order transitions separates phases 
of different density. For the continuous spin ice model we have mapped out the equivalent 
phase diagram by Monte Carlo simulation. In each case the phase boundary was estimated 
from the maximum in the magnetic susceptibility. In principle, now we know the true ground 
state, we could measure the true order parameter and its fluctuations, by calculating the 
projection onto the calculated state. 

The results are shown in Fig. |]a, and represented schematically in Fig. ||b. It is seen that 
the two first-order lines of the spin ice phase diagram coalesce below the ordering temperature, 
which is therefore seen to be a triple point (see Fig. Hb). The triple point temperature increases 
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Fig. 3 - Comparison of the numerical zero temperature magnetization mo, estimated by a zero- 
temperature Monte Carlo method, with that estimated by the analytical expressions ([]) (left hand 
line) and (^|) (right hand line). Both asymptotic expressions are seen to break down near the tricritical 
point D t J J = 5.0 ± 0.5. 



with decreasing D / J and eventually the two "wings" disappear and give way to a line of first- 
order transitions along the zero field axis, typical of a ferromagnet. We would anticipate that 
the two wings disappear precisely at the tricritical point determined above, D tc / J ~ 5. This 
kind of phase diagram is exhibited by idealised ferroelectrics and also some real ones such 
as BaTiOa It can be rationalised by Landau theory |23|, ^JJ, where the free energy is 

written G = Go + ^-m 2 + ^fm 4 + ^m 6 — NmH. Here, m and H are the order parameter and 
conjugate applied field and the c„'s are constants related to the n-th order susceptibilities. If 
C4 > then the zero-field transition is second-order and if C4 < it is first-order, with the 
characteristic "winged" phase diagram. The slopes of the wings are approximately ~ \J~%- I n 

general, termination points of first-order lines need not be associated with critical fluctuations, 
in which case they should be termed "first-order" critical points. One such point is known 
to occur on the zero-field axis of the phase diagram of the one dimensional Ising ferromagnet 
with additional inverse square interaction ^6|. In fact it has been suggested j24| that the 
end-point of a line of symmetry sustaining transitions is usually first-order, the liquid gas 
transition being a special case. In general, at a first-order critical point, the susceptibility will 
diverge, but there may or may not be critical fluctuations ^4j,^5|). A detailed analysis of this 
question with regard to the current system would be an interesting topic of future study. 

Conclusions. - In conclusion, it is noteworthy that the continuous spin ice model displays 
behaviour characteristic of ferroelectrics. It is plausible that ice itself would exhibit such a 
phase diagram if the dynamics did not become immeasurably slow at temperatures well above 
the triple point ||| . However ice is not proton ordered in the absence of defects or stabilizing 
surfaces and there is much debate concerning the true nature of its ground state p7j-|3l[|. 
Very similar remarks apply to Ho2Ti20y, where the dynamics also become immeasurably 
slow below 0.7 K. In view of our results, it might be that the ordered state favoured by the 
long range part of the dipolar coupling competes with that favoured by continuous spins. 
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Fig. 4 - Magnetic phase diagram for the continuous spin ice model define by eqn. (|l|). (a) Phase 
diagram determined by Monte Carlo simulation on a system of N — 2000 spins. Points on the coex- 
istence lines were determined from the maximum in the susceptibility in fixed temperature (vertical 
error bars) or fixed field (horizontal error bars) Monte Carlo runs. For ease of comparison with ex- 
periment and ref.( |H|), the field h is defined as h = H x (/^b/&b = if/0.06717, with /ib and fee the 
Bohr magneton and Boltzmann's constant respectively, and H defined in eqnQ (b) Schematic phase 
diagram showing the case of D/J = (short dashed line), D/J = Dtc/J (thin line), D/J > Dt c /J 
(thick line) and D/J — > oo (long dashed line) 



Through the application of a magnetic field in the [1,1,0] direction S, it is possible to put 
Ho2Ti207 into an ordered state of the same symmetry as the ground state of the dipolar spin 
ice model pij ] . It would be interesting to prepare a sample in this way at a temperature below 
the expected ordering transition at 0.18 K, and to examine the stability of the ordered state 
as the field is removed. At a theoretical level, it would be extremely interesting to consider 
the effect of continuous spins on the dipolar spin ice model of den Hertog and Gingras p5[ , 
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that is an accurate description of Ho2Ti 2 07 and Dy 2 Ti 2 07 |l5|,|l8|]. 
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